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Abstract. We consider the following eigenvalue optimization problem: 
Given a bounded domain SI C R" and numbers a > 0, A £ [0, find 
a subset D C of area A for which the first Dirichlet eigenvalue of the 
operator —A + ol\d is as small as possible. 

We prove existence of solutions and investigate their qualitative prop- 
erties. For example, we show that for some symmetric domains (thin 
annuli and dumbbells with narrow handle) optimal solutions must pos- 
sess fewer symmetries than Q; on the other hand, for convex fl reflection 
symmetries are preserved. 

Also, we present numerical results and formulate some conjectures 
suggested by them. 



1. Problem and Main Results 

We study qualitative properties of solutions of a certain eigenvalue opti- 
mization problem. In physical terms, the problem can be stated as follows: 

Problem (P) Build a body of prescribed shape out of given 
materials (of varying densities) in such a way that the body has a 
prescribed mass and so that the basic frequency of the resulting 
membrane (with fixed boundary) is as small as possible. 

In fact, we will consider a more general problem, which we now state in 
mathematical terms: Given a domain $7 C M. n (bounded, connected, with 
Lipschitz boundary) and numbers a > 0, A G [0, (with | • | denoting 
volume). For any measurable subset D C £7 let \D be its characteristic 
function and X^i(a,D) the lowest eigenvalue A of the problem 

—An + avnii = Xu on 17 

(1) 

K J u = o on an. 

Define 

(2) A n (a, A) = krf An(a,D). 

\D\=A 
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Any minimizer D in (j2|) will be called an optimal configuration for the data 
(Q,a,A). If D is an optimal configuration and u satisfies (p]) then (u, D) 
will be called an optimal pair (or solution). Our problem now reads: 

Problem (M) Study existence, uniqueness and qualitative prop- 
erties of optimal pairs. 

As is well-known, u is uniquely determined, up to a scalar multiple, by D, 
and may be chosen to be positive on Q. In addition, we will always assume 

u 2 = l. 

n 

(Integrals over SI are always taken with respect to the standard measure.) 
Clearly, changing D by a set of measure zero does not affect Ao(a, D) or u. 
Therefore, we will consider sets D that differ by a null-set as equal. 

At first sight, it is not obvious that problem (M) generalizes problem (P). 



In fact, we will see (Theorem 13) that there is a number a^(A) > such that 
solutions of problem (P) are in one to one correspondence with solutions of 
problem (M) with parameters in the range a < ao,{A). The number an(A) 
is characterized as the unique value of a satisfying 

(3) A n (a n (A),A) = a n (A), 



see Proposition 10 



Our investigations are theoretical and numerical: Numerical results (ob- 
tained by M.I. and I.O.) suggest properties of optimal configurations; this 
leads to the formulation of conjectures, and some of these are proved rigor- 
ously (by S.C., D.G. and K.K.). 

A central tool in our investigations is the variational characterization of 
the eigenvalue: 

fo |Vn| 2 + a L ydu 2 
X n {a,D)= inf Rn(u,a,D), Rn(u,a,D) - " v 



uefl-i(n) J s 



and the eigenfunction u is a minimizer. So An(a,A) is characterized by 
An(a,A)= inf Rn(u,a,D). 

\D\=A 

We first prove the following theorem on existence and basic properties of 
solutions. It is fundamental for all further considerations. 

Theorem 1. For any a > and A € [0, |fi|] there exists an optimal pair. 
Moreover, any optimal pair (u, D) has the following properties: 

(a) u G C l ' 5 (Q) n H 2 (n) n C 7 (H) for some 7 > and every 6 < 1. 

(b) D is a sublevel set of u, i.e. there is a number t > such that 

D = {u< t}. 

(c) Every level set {u = s}, s > 0, has measure zero, except possibly in the 
case a = an(A), s = t. 
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Here we use the short notation {u = t} = {x : u{x) = t}. Since xd is 
discontinuous, solutions u may not be twice differentiable, so equation (||) 
is understood in the weak sense. 

Note that Theorem |l|(b) shows in particular that our problem is equiva- 
lent to finding the smallest eigenvalue and associated eigenfunctions of the 
nonlinear problem (with free variables u and t) 

—Au + ax{ u <t} u = ^ u on ^ 
(4) u = on dn 

\{u < t}\ = A. 

The question of uniqueness is much more subtle: For some domains f2 
there will be a unique optimal pair for all a, A, while for others there will be 
many, for certain ranges of a, A. This follows from our results on symmetry 
preservation and symmetry breaking below. 

We now list a few questions that naturally come to mind: 

(SY) If f2 has symmetries, does D have the same symmetries? (Note that 
if 0, and D have a symmetry in common then u will also have this 
symmetry since it is uniquely determined by f2 and D.) 

(CX) Assume £1 is convex. Is D c := f2 \ D convex? Is D unique? 

(CN) Is D or D c connected? 

(FB) What is the regularity of the free boundary dD? 

We give partial answers to all of these questions. Some proofs, mainly 
relating to (FB), and additional results can be found in the companion paper 
[CGK]. Many open problems remain, see Section^. 

At this point, the reader is invited to look at Figures 1-3 for a first im- 
pression. 

We now state our qualitative results. As a general convention, constants 
only depend on the quantities indicated as subscripts or in parentheses, 
unless otherwise specified. Often we suppress the subscript Q. 

First, as an easy consequence of Theorem |l| one has: 

Theorem 2. Fix a > 0, A > 0, and let D be an optimal configuration. 

(a) D contains a tubular neighborhood of the boundary dQ. 

(b) If a < afi(A) then every connected component Dq of the interior of D 
hits the boundary, i.e. Dq n <9£1 ^ 0. 

In particular, if O is simply connected and a < an(A) then D is connected. 

The number an (A) was defined above, see (0). The significance of the 
condition a < a$j (A) is that it is equivalent to A-u < on Q. One always 
has 

an (A) > fin. 

Here and throughout the paper, fin denotes the first eigenvalue of the Dirich- 
let Laplacian on Q, and ipn the positive, L 2 -normalized eigenf unction: 

-A*l>n = Unipn on $7, ipn = on d£l. 
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Next, we consider the dependence of Aq and solutions (u, D) on a and A. 
Here it is convenient to formulate our problem also for a = 0, as follows: If 
a = then a solution (unique in this case) is a pair (ijjfi,D) where D is the 
sublevel set of i/jq of area A. (Since i/jq is real analytic and non-constant, 
such D exists for every A and is unique.) 

We will prove strict monotonicity and Lipschitz continuity of in both 
parameters (Proposition |l0|). Continuous dependence of optimal pairs (u, D) 
on the parameters may be expected only at parameter values where they are 
unique. This is the case, in particular, if a = or A = or A = in these 
cases u = ipn, and the continuity is proved in [CGK]. Here we only state 
the results. They are used only in the proof of Theorem ||. For example, we 
have the following: 

Theorem 3. For s > let [Q] s = {ipn < s}, where tpn is the positive L 2 - 
normalized first eigenfunction of —A on f2. Fix A £ [0, |0|] and choose tci 
such that \[£i\ tn \ = A. Then for any 6 > there is «o = &o(8,Q) such 
that whenever a < «o an d D is an optimal configuration for (a, ^4) then 
\t — tn\ < 5 and 

CD C [n] tn+S . 

We now address questions of symmetry. First, we prove symmetry preser- 
vation in the presence of convexity: 

Theorem 4. Assume that the domain £1 is symmetric and convex with re- 
spect to the hyperplane {x\ =0}. In other words, for each x' = (x%, . . . , x n ) 
the set 

(5) {xi : (x 1 ,x') € ft} 

is either empty or an interval of the form (— c, c). 

Then for any solution (u, D) both u and D are symmetric with respect to 
{x\ = 0}, D c is convex with respect to {x\ = 0} ; and u is decreasing in x\ 
for x\ > 0. 

For example, any solution in an elliptic region has a double reflection 
symmetry, see Figure 1. The principal tool here is Steiner symmetrization. 
See [K2] for an overview on such methods. Theorem || easily implies the 
following uniqueness result (the only case where we can prove uniqueness!): 

Corollary 5. Let fi = {|x|<l}&e the ball. Then there is a unique optimal 
configuration D for any a, A, and D is a shell region 

D = {x : r(A) < \x\ < 1}. 

One of the most interesting phenomena studied in this paper is symmetry 
breaking for certain plane domains f2. That is, an optimal configuration D 
may have less symmetry than fi. We will prove it for two types of domains: 
Thin annuli and dumbbells with narrow handle. An annulus has rotational 
symmetry, a dumbbell has a reflection symmetry. 
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Theorem 6. Fix a > and 5 G (0, 1). For a > let Q a = {x G K 2 : a < 

|x| < a + 1}. There exists ao = ao(a,<5) suc/i i/toi whenever a > ao and D 
is an optimal configuration for £l a with parameters a and A = o~|f2 a | £/ien D 
is not rotationally symmetric. 

See Figure 2. For dumbbells we prove a little more than symmetry break- 
ing: 

Theorem 7. For h G (0, 1) define the dumbbell with handle width 2h 

fl h = B 1 (-2, 0) U ((-2, 2) x (-h, h)) U B 1 (2, 0) 

where B r (p) = {x G M? : \x — p\ < r}. Fix a > and A G (0, 2tt). Then 
there is h$ = ho(a, A) > such that we have for h < h^: 

(a) Any optimal pair (u,D) is not symmetric with respect to the x^-axis. 

(b) If A > tt then for any optimal pair (u, D) the complement D c is con- 
tained in one of the lobes (i.e. one of the balls -E>i(±2, 0)). 

See Figure 3. In fact, similar results hold for more general dumbbells. 

As we remarked before, symmetry breaking implies non- uniqueness: For 
example for a dumbbell the pair («', D') obtained from a solution (u, D) by 
reflection in the X2-axis will be a solution, and different from (it, D) by the 
theorem. 

The following result on the regularity of the free boundary is proved in 
[CGK]: 

Theorem 8. If (u, D) is an optimal pair, x G dD and Vn(x) ^ then dD 
is a real analytic hypersurface near x. 

The difficulty is that xd is discontinuous at x G dD, so u is not even C 2 
there. That the level set {u = t} has C u regularity nevertheless is proved 
by introduction of suitable local coordinates (with u as one coordinate) and 
analysis of the resulting nonlinear elliptic equation. 

Similar arguments and continuity considerations for a near zero allow us 
to give partial answers to problems (CX) and (FB): 

Theorem 9. Suppose Q is convex and has a C 2 boundary. Then there is 
ao(A, Q) > such that for any a < ao an d an U optimal configuration D, 
one has: 

(a) dD flH is real analytic. 

(b) D c is convex. 

Problem (P) and generalizations of it (to higher eigenvalues and to a max- 
imization problem), but with fewer qualitative results, were studied before in 
[Kr], [CM], and [C] (where Theorem |] is stated, but the proof is incomplete 
since the case of equality in the rearrangement inequalities is not addressed) . 

Problems similar to problem (M) (e.g. with IP potentials) were considered 
in [AH], [Eg], [AHS], [CL], and [HKK] . 
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The paper is organized as follows: In Section @, we prove Theorems |] 
and H and discuss the parameter dependence of Aq. Also, in Subsection 
2.3 we discuss the relation of problems (P) and (M). In Section |3| we prove 
Theorems ^, and [7 on symmetry questions, and Corollary |!| In Section 
|] we prove Theorem 9|. In Section |5| we describe the numerical algorithm 
used. In Section |6] we state some open problems and conjectures. Finally, 
we collect some standard facts about elliptic PDEs in the Appendix. 
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2. Basic results 

2.1. Existence and regularity. Proof of Theorem El We first prove 
existence and regularity: The regularity statements in (a) hold for solutions 
of equations 

— Au + pu = 

with p bounded by standard elliptic theory, see for example [GT, Theorem 
8.29 and Corollary 8.36]. 

To prove existence, fix a and A, and write A = A^i(a,A), \(D) = 
Aq(ci,D) for simplicity. Let Dj be a minimizing sequence, i.e. X(Dj) —* A 
as j — ► oo. Let Uj € Hq (all function spaces are defined on Q) be the 
positive L 2 -normalized first eigenfunction of —A + ctXDy Since \(Dj) is 
bounded, the sequence {uj} is bounded in Hq. Also, {xDj} is bounded in 
L 2 . Therefore, we may choose a subsequence (again denoted Uj,Dj) and 
u € Hq, r] € I? such that uj — 1 u in Hq (weak convergence) and xd 3 V i n 
I? . This implies Uj — ► u (strongly) in L 2 , XD 3 Uj — 1 V u m -^ 2 > an d In 7 ! = A- 
Now taking limits in the weak form of the eigenvalue equation 

/ Vuj ■ Vtp + a XDjUjip = A(-Dj) / Ujtp Vip £ Hi 
Jn Jn Jn 

we get 



(6) 



—Am + ar/u = Au (weakly). 



SYMMETRY BREAKING FOR COMPOSITE MEMBRANES 



7 



We have 

< r] < 1 a.e. 

since < xd < 1 f° r all j and weak convergence preserves pointwise in- 
equalities a.e. (exercise!). Therefore, u has the regularity stated in (a). 

It remains to prove that r\ may be replaced by a characteristic function. 
Since f n u 2 = 1, (||) shows that 

(7) I \Vu\ 2 + a [ r/u 2 = A. 

Jn Jn 

Now the minimization problem 

inf / rjv 2 
vj v=AJn 

0<7/<l 

has a solution r/ = \D where D is any set with \D\ = A and 

(8) {u < t} C D C {u < t}, t : = sup{s : \{u < s}\ < A} 

(compare the 'bathtub principle', Theorem 1.18 in [LL]). Therefore, we get 
from (0) 

/ |Vn[ 2 + a xdu 2 < A. 
Jn Jn 

By definition of A as a minimum, this must actually be an equality, and 
(u, D) is a solution. 

(b) Let (u, D) be any solution. Then it is obvious that (||) must hold 
(always up to a set of measure zero; if @) didn't hold then one could reduce 
J D u 2 by shifting a part of D from {u > t} to {u < t}). Set M s = {u = s} 
for any s > 0. Using Lemma 7.7 from [GT] twice, we see that Au = a.e. 
on J\f s (since u = const on J\f s ; recall that u is in H 2 ). Therefore, 

(9) (A — oxd)u = a.e. on jV s . 

Since u > and A > 0, this shows that D° n M s has measure zero. Taking 
s = t we get (b). 

(c) If s > t then 7V S C D c , so \M S \ = by (^). The same argument works 
if s = t and q/A. 

Finally, u satisfies — Au = (A — a)u on the open set {u < t}, hence u 
is real analytic there, and therefore the level sets Af s have measure zero for 
s <t. □ 



Proof of Theorem |]. Part (a) is clear from Theorem |Tj(t>) . To prove 
(b), assume this was false. Then there is an open subset Dq C {u < t} 
with 8Dq C D c = {u > t} and therefore u = t on ODq. Then u assumes 
a minimum at some xq £ -Do- But this is a contradiction since a < an(A) 
implies A(a, A) > a (see Proposition |l0| below) and therefore Au = (a — 
A(a, A))u < on D . □ 
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2.2. Parameter dependence of A. 

Proposition 10. (a) The function (a, A) \— > A(a, A) is Lipschitz contin- 
uous, uniformly on bounded sets. More precisely, we have, for any 
a, a' > 0, A, A' € [0, 

(10) \A(a,A) - A(a',A')\ < 

l« - a I |^j A \ mmja, a }C n ^ {oiya , } 

with Co, Q bounded for a bounded. 

(b) A(a,A) is strictly increasing in A for fixed a > 0, strictly increasing 
in a for fixed A > 0, and A(a, A) — a is strictly decreasing in a for 
fixed A < |0|. 

(c) If A < |f2| i/ien i/iere is a unique value a = oq(A) with 

(11) A(a n (^),A) = a n (A). 

The function 7Sq is continuous and strictly increasing, an(0) = //n, and 
c?n(A) — > oo as ^4 — > 



Proo/. (a) Write A = A(a, A) and A' = A(a', A'), and let (u, D), (u' , D') 
be minimizers for A, A' respectively. We may assume JqU 2 = Jq(u') 2 = 1, 
so that 



A= f |Vu| 2 + a f u 2 , \D\ =A, 

JCl JD 



and similarly for A' etc. By symmetry of ([!(]) we may assume that A' > A. 
Choose D\ C D' with \D\ \ = A and D[ D D with \D[\ = A'. Here we may 
assume that D[ is of the form {u < s} for a suitable number s. Using the 
optimality of (u, D) for A we get 

(12) A< f \Vu'\ 2 + a [ (u') 2 = A' + (a - a') [ (u') 2 - a [ (u') 2 . 

Jfl JD 1 JD 1 J D'\D 1 

Similarly, using the optimality of (u' , D') for A' we get 

(13) A'< \Vu\ 2 + a u 2 = A + {a'-a) u 2 + a 

Jn Jd[ Jd' l Jd'Ad 

Alternatively, we may rewrite this as 



u 2 . 



(|f) A' < A + (a - a) \ u 1 + a' I u 



2 • n' / 2 

D Jd[\D 



In order to estimate the integrals in (|12D, (|13|) and (|13r) which are multiplied 
by ±(a — a'), observe that for any s > and any function u we have 

ku<s} u2 \{U<S}\ 

f n u 2 - \n\ ■ 
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The other integrals are estimated using the uniform estimate ( f47| ) : u solves 
the equation — Au + axD u = An. A is bounded in terms of f2 and a since 
one may apply ( |i~2| ) with a' = 0, A = A', to obtain A < + a. Therefore, 
the uniform bound (|47|), applied to G = Q, yields 

[ u 2 < (A' - A) sup u 2 < (A' - A)Cn, a . 
Jd[\d n 

Finally, we obtain ( |Tc|) by applying these estimates to (|l2|) and (jl3|) in the 
case a < a', and to (|12|) and (|l3|') if a > a'. 

(b) This follows immediately from (|l^) and the unique continuation the- 
orem. 

(c) This follows easily from (a) and (b) since A(a, A) — a equals fin > 
for a = and tends to — oo as a — > oo by (a). □ 

We now consider continuous dependence of optimal pairs (u, D) on the 
data. First, near a = 0: 

Proposition 11. Fix D C £1. Let u a ^r> be the (positive, L 2 -normalized) 
first eigenfunction of — A + a\D, and ipa = no,D the first eigenf unction of 
—A. Then there is a constant C = Cn such that, for < a < 1, 

\\u a ,D < Ca, 

in the H 2 (n) and L°°(0) norms, and in C^ S (Q) if dU is in C 1 ' 5 ^). 



Proof See [CGK]. □ 

Proof of Theorem [|. This is almost immediate from Proposition [ll|, see 
[CGK]. □ 
Similarly, one has continuity in A at A = and at A = Here we only 
consider the latter case: 

Proposition 12. Let Q be a smooth bounded domain and fix a > 0. Let 

M = maxipQ. 

Then, for any 5 > there is Aq = Aq(6, a, O) < |0| such that whenever 
A > Aq and D is an optimal configuration for (a, A) then 

D c c {Vn > M — 5}. 



Proof See [CGK]. □ 

2.3. Relation of problems (P) and (M). We want to show that problem 
(P) (see Section [j]) is a special case of problem (M). 
The mathematical formulation of problem (P) is: 

Given < h < H (lower and upper bounds for the densities of the materi- 
als that are available) and the prescribed total mass M £ [h\tt\, -ff|0|], M > 
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0, consider measurable 'density functions' p satisfying 

h< p<H, p = M. 
Jn 

Then the objective is to find p and u which realize the minimum in 

fa |Vn| 2 

(14) e(h,H,M) :=inf inf 

P u£HZ(n) J n pu 2 

The corresponding eigenvalue problem is 

(15) -An = Qpu, u\ dn = 0. 

(We assume the modulus of elasticity to be the same for all materials.) 
Problem (P) and problem (M) are related in the following way: 

Theorem 13. (a) If (u,p) is a minimizer for problem (P) then p is of the 
form 

PD = hxD + HXD? 
for a set D of the form D = {u < t}. That is, only two types of 
materials occur. 

(b) The pair (u, po) is a minimizer for problem (P), with parameter val- 
ues (h,H,M), if and only if (u,D) is a minimizer (optimal pair) for 
problem (M), with parameter values (a, A) given by 

(16) a = (H — h)@(h, H, M), 

(17) A = ^ihr- 

The minimal eigenvalues are related by 

(18) A(a, A) = H@(h, H, M). 

(c) The values of (a, A) that occur when h, H, M vary are precisely those 
satisfying 

Ae[0, \n\), 0<a<a n (A) or 
A=\Q\, < a < co, 



where an(A) is defined in (|il|). In particular, a = an(A) corresponds 
to h = 0. 

Note that problem (P) really depends on two parameters only since for 
k > one has 

B(nh, kH, kM) = /c _1 e(/i, H, M), 
with the same minimizers (up to a factor k for p). This is obvious from (14). 

Proof. (a) This is almost obvious from (|i~4]), and proved just like part (b) 
of Theorem |l[ 

(b) First, if p = p D and \D\ = A then M = J Q p = Ah + (|0| - A)H, 
which gives (|l~7|). 
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Simple manipulation shows that 

(19) -An = Qp D u = 0(h XD + Hxdc)u 
is equivalent to 

(20) -Au+(H -h)Q X DU = HQu. 

Now if (u, pr>) is a minimizer for problem (P) then it satisfies ( |l9|) with 
G = e{h,H,M), and then @ shows that A(a,A) < HO{h,H,M) with a 
satisfying (|l6|). 

Conversely, if (u, D) is a minimizer for problem (M) with parameter values 
(a, A) given byjl|), (|l7|) then @ holds with HQ replaced by A = A(a, A), 
so instead of (|9) we get -Au = Qp D u+(A-HQ)u where 6 = Q(h, H, M). 
Multiplying by u and integrating gives 

[ |V-u| 2 = 9 f p D u 2 + (A - HQ) [ u 2 . 
Jn Jn Jn 

Now the definition of Q implies that \X7u\ 2 > Q J^pdu 2 , so we get A > 
HQ. This proves A(a, A) = HQ(h,H,M) and part (b). 

(c) If A = then D = O, p = h and therefore hQ(h, H, M) = p,^ from 
(|l5|) , so a = pn can take any positive value by suitable choice of h and 
H. 



Now let A < |f2[. By Proposition 10(b) and (c), a varies in the indicated 



range precisely when A(a, A) — a varies in 

[0,pn)- From (||) and @ one 
has 

A(a, A) - a = := /i6(/i, H, M), 

so we only need to show that hQ has range [0,pn) (with A fixed). First, 
hQ > by definition, and hQ = A — a < pn by Proposition |l0j, since 
a = (H — h)Q > 0, so the range of hQ is contained in [0, pn)- Next, hQ = 
for h = (and then M can be adjusted to ^4), and in the limit H = h one 
has p = h and hQ = pn, so when H — > /i then /iB — > u^, and clearly M 
can be adjusted to A. Using continuity of hQ (which is proved as for A in 
Proposition [Toj) we get the claim. □ 



3. Symmetry preservation and symmetry breaking 

3.1. Symmetry preservation in the presence of convexity. Here we 
prove Theorem |J 

Proof of Theorem [J. We use Steiner symmetrization (symmetrically 
decreasing rearrangement) u i— > u& with respect to the hyperplane {x\ = 0}. 
This is defined as follows. Assume u £ Hq(£1) nC°(J7): For each x', u#(-, x') 
is the unique function of xi which is symmetric in x\ and decreasing for 
x\ > such that |{xi : u^{x\,x') > t}\ = \{xi : u{x\,x') > t}\ for all t £ R. 
It is well-known (see, e.g., [LL], [AB]) that, for all x' and i = 1,... with 
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integrals taken over the set (|), 

(21) J \d Xi u#\ 2 dx x < J\d Xt u\ 2 d Xl , 

(22) J {u*f dx x = J 



u dx\, 



(23) J (o>xd)#(u*) dxi < J axDU dx x . 

Here, is the increasing symmetric rearrangement of a function /, which 
is defined by = — (—/)*. Note that ( ^T| ) for i = 1 is just the standard 
rearrangement inequality in one dimension, while for i > 1 it is proved as 
follows: Replace the partial derivatives by difference quotients (v e (xi) — 
vo(xi))/e with v e {x\) = u(xi, ... , Xi + e, ...). After multplication by e 2 the 
claimed inequality becomes simply J \vf — v^\ 2 dx\ < J \v e — vo\ 2 dx\ which 
is well-known. 

Fix q and A and assume (u, D) is an optimal pair. Define the set by 
Xd* = (xd)#- Integrating ([H]), (^) and ( |23| ) over x' and summing (21) 
over i we get 

# J n |Vn#| 2 Q ;x + / n (axD)#(n # ) 2 dx 

(24) " J^?& -A(a,D). 

Since we have \D&\ = \D\ = A (by (p2^) applied to Xd), optimality of (u, Z?) 
implies that (vft, D&) is also a minimizer and that equality holds in ( ^Tl) and 
(^3|), for all i and almost all x' . We need to show that this implies u = it*. 
The statements about D then follow from the characterization D = {u < t}. 

First note that since (u#,D#) is a minimizer, the function u# solves 
the equation — + axn#u^ = \(a;D#)u#. Therefore, u and are 
continuously differentiable by Theorem [l], so equality in (21) holds for all 



x 1 . By a result of Brothers and Ziemer (see [BZ]) this equality implies 
vft{x\,x') = u(xi,x') for all x\ provided the set \x\ : d Xl vft(xi,x') = 0} 
has measure zero. 

Therefore, we will be done once we have shown that the set 

(*) {v = 0} has measure zero, where v = d Xl u^. 

We will give two proofs of this: The first proof works whenever a ^ an(A) 
and the second proof works whenever a < an (A), so together they cover all 
cases. 

First proof of (*), assuming a ^ olq,(A): Assume this was not so. Define 
i# by = {ti* < t&}. v satisfies — Av + axn# v = \{a,D#)v on {u# / 
i*}. Since {u^ = t^} has measure zero by Theorem |l| and the assumption 
a 7^ &n{A), v vanishes on a set of positive measure in the open set {u* ^ 
ffcy, so the unique continuation theorem (for sets of positive measure, see 
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[FG]) applied to v implies that v = on some connected component K 
of 7^ t*}. Therefore, vft is constant in the xi-direction on K. Since 
v# = or on dK we conclude that then vft must actually be constant 
on K. This is a contradiction to Theorem [l|(c). 

Second proof of (*), assuming a < a^i(A) (this proof is taken from Cox 
[C]): We show that actually v < for x\ > 0, so that {v = 0} is contained 
in the hyperplane {x\ = 0}. We have — Au* = A(a,A)u# — axD# u ^i 
and the right hand side is decreasing in x\ (for x\ > 0) by definition of 
the rearrangement and since a < A(a, A) by Proposition Taking the 
xi-derivative (in the sense of distributions), we get Av > as distribution. 
Also, v is continuous, so by the classical theory of subharmonic functions 
it satisfies the maximum principle (alternatively, it is in H 1 and then the 
maximum principle as in [GT], Ch. 8, applies). Since v < 0, we conclude 
that v < unless v vanishes identically in x\ > 0, which is clearly impossible. 
This proves (*). 

This concludes the proof that u = u& and hence the proof of the theorem. 
Note that in the case a < an(A) the second proof of (*) above actually shows 
that u xi < for x\ > 0. 

□ 

Proof of Corollary The only set D C {\x\ < 1} which has the 

symmetry and convexity properties stated in Theorem |] in all directions is a 
shell region as stated. Clearly, r(A) is uniquely determined by A. Therefore, 
D is unique. □ 

3.2. Symmetry breaking on annuli. We now give the proof of Theorem 
H about symmetry breaking on an annulus, 

= Q a = {x E R 2 ; a < \x\ < a + 1}, a > 0. 

Let D be any radial set in Q, 

D = {(r,9);r G D u < 9 < 2vr}, D x C (a,a + l), 

and let u be the first eigenfunction for D, with eigenvalue a: 

(25) —Au + axDU = au onQ, u\qq = 0. 

For a sufficiently large (depending on a and 5 = \D\/\£l\) we will construct 
a comparison domain D and a function u which satisfy 



In a l v ^l 2 + Iu a Xd 



u 2 I 



(26) - r " a < a. 

This shows that D is not an optimal configuration and hence implies the 
theorem. 

In order to construct D and u, first pick N = N(5) with 

1 



14 S. CHANILLO, D. GRIESER, M. IMAI, K. KURATA, AND I. OHNISHI 

and consider the sector 

E + = n a n {(r, 0); < < n/N}. 

Then let u be the first Dirichlet eigenfunction of the Laplacian on E + and 
\\(E + ) be the first eigenvalue, 

(27) -Au = Xi(E + )u on E+, u\ 9E+ = 0, 

extended by zero on Q\E + ; the set D can be taken to be any subset of £l\E + 
with \D\ = \D\. This is possible since \D\/\Q\ = 6 < I - ^ = \Q \ E + \/\Q\. 

Note that since suppuHD = 0, we have (J n |Vu| 2 + f n Xd^ 2 )/ f n u 2 = 
Ie l^ 7 ^! 2 / Ie ^ 2 = so (^) * s equivalent to 

(28) Ai(£+) < a. 

In order to prove this, we need to introduce a third eigenvalue problem, 
which is intermediate between ( p5| ) and (p7|). 

Define v to be the lowest eigenfunction for the problem ( p5|) among func- 
tions of the form 

v(r,6) = h(r) sin NO, 

and let r be the associated eigenvalue. Note that problem ( |25| ) for such 
functions is equivalent to the problem 

1 N 2 

(29) -h"(r)--ti(r) + - T h(r) + axD 1 (r)h(r) = rh(r) on[a,a + l], 

(30) h{a) = h(a + 1) = 

for h. Thus, h is the first eigenfunction of this Sturm-Liouville problem, and 
the eigenvalue t is characterized by 

(31) t= inf Ja VW ' 1 '- ' 



where S = {g G C^a, a + l];g(a) = g(a + 1) = 0}. From this the (well- 
known) fact that h does not change sign on [a, a + 1] is evident; so we may 
assume 

h> 0. 

We will compare u with v and t> with u. The following two lemmas provide 
the needed estimates. 



Lemma 14. Let a be the lowest eigenvalue for the problem (25) (with D 
radial) on fi a & = {x € M? : a < \x\ < b}, and let r be the lowest eigenvalue 
for eigen] 'unctions of the form v(r, 9) = h(r) sin N9 on fi a 6- Then we have 

r-a< N 2 /a 2 . 
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Proof. Since xd is assumed radial, the first eigenfunction of ( |25| ) is a radial 
function u = f(r). Now consider the trial function w(r,9) = f(r)smN6. 
We have 

fn a ^ Vw ^ + aXDW 2 )dx 
f n w 2 dx 

Thus, 

Ia((f(r)) 2 + ^f(r) 2 + a X Dj(r) 2 )rdr 



T < 



By definition of f(r) we get 



/a f( r ) 2r dr 



J a f(r) 2 rdr 

The claim follows. □ 

Lemma 15. Define v as above. Assume D is radial and \D\/\Q\ = 5. There 
exists a positive constant c a ^, independent of a, such that for all a > 1 we 
have 

fr, V 2 dx 

±il > c x 

Jn v dx 



Proof. We see from v(r,8) = h(r) sin N9 that 

f D v 2 dx = J a a+1 xp^Krfrdr 
L v2dx C'Hrfrdr 

h satisfies equation (|29|). For r one has a uniform bound r < C a ^ with C a j 
independent of a > 1, because from ( j3l"l ) one gets 

f a+1 (g') 2 rdr 

t < inf Ja a+ ^ 1 + a + N , 

J a g 2 r dr 

and by using for g the translate of any fixed test function on [0, 1] one sees 
that the first term on the right is bounded by some absolute constant. 

Therefore, the coefficients of equation (^) are uniformly bounded for 
a > 1. Also, we have h > 0. Lemma 16 in Section |6| then implies that one 
has 

(33) LX/ M, xl h>c a ^\\h\\ L 2 { +1) . 

[a+6/4,a+l~5/4] v ' 



Since \Di \ = 5, we have \[a + 5/4, a + 1 - 5/4} f\D x \ > 5/2. Therefore, 
ra+l $ 

(34) / XDx{ r )h{ r ) 2r dr > —a inf h 2 

Ja 2 [a+S/4,a+l-8/4] 

and 

ra+l ra+l r 

(35) / h(r) 2 r dr < (a + 1) h 2 <2a 

J a J a J a 



h 2 . 
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Combining d33|), (j34[) and ( ^5| ) with (^) we get the Lemma. □ 

End of proof of Theorem We have 

r = Jnl^l 2 ^ + aJ n XDV 2 dx 
fn v2dx L v2dx 
Since v(r, 6) = h(r) sin NO, v vanishes on the rays 8 = and 9 = ir/N. Since 
\v\ and |Vu| are periodic in 6 of period n/N, we can replace f2 by E + in 
the first quotient. Therefore, we can use v as test function in the Rayleigh 
quotient for the Dirichlet Laplacian on E + and obtain 

L\Vv\ 2 dx J E+ \Vv\ 2 dx ^ 
J n v 2 dx f E v 2 dx 

Combining this with (|36| ) and Lemma [l5| we therefore get 
(37) r > Ai(^+) +ac a , 5 . 



From Lemma |lj we then get 

cr > r - N 2 /a 2 > Xi(E + ) + ac a>5 - N 2 /a 2 . 

If a is chosen so large that N 2 /a 2 < ac a ^ then this gives (|28| ) and hence the 
theorem. □ 



3.3. Symmetry breaking on dumbbells. Proof of Theorem [?|. Since 
a is fixed throughout, we will write Xn(D) = Xn(a, D), Aq(A) = A^(a, A) = 
inf|£>| = ^ Xq(D). Here we keep the index Q since we will also consider these 
quantities with £1 replaced by one of the 'lobes' B± = i?i(±2,0). All (im- 
plied) constants will only depend on a and A. Write A B = A B± , and given 
D, let 

D± = Dr\B±, A± = \D±\. 

Further, we introduce 

A min = min{min(|L>_|, \D+\) : D C n, \D\ = A}. 

Thus, if D is distributed over f2 with the greatest possible imbalance between 
D + and D_ then the smaller of D± will have area ^4 m in- It is easily checked 
that 

^4min = max(0, A - \B C _\). 
We first sketch the idea of the proof: 
1. For h = 0, i.e. two disconnected balls, one clearly has 

(38) A n (A) = mm(A B (A^),A B (A + )). 

Since A^ is strictly increasing, it is optimal to put as much of D as possible 
in one ball, say B + , and the 'small' remainder in the other. Thus 

A n (A) = A B (A min ), 

and the eigenfunction is zero in B + . 
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2. For small positive h, this situation should be approximately the same: 
Equation ( |38| ) will hold with an error that is a power of h (compare equation 
( f42| ) below), so the same argument as in 1. implies symmetry breaking. Also, 
the eigenfunction must be small on one lobe, and since D = {u < t}, one 
gets (b) from an estimate of t. 

We now carry out the details. Let (u, D) be an optimal pair. Assume 

u 2 = l. 

n 

First we need an estimate ensuring that the perturbation introduced by 
the handle is small. This is provided by the following estimate near the 
boundary (see [GT, Theorem 8.27 with Rq = 1 and R = 3ft]), which is 
applicable since O satisfies a uniform exterior cone condition (uniformly in 
ft): There is (3 G (0, 1] such that 

(39) max u(x) < Ch^\\u\\ L 2 (Q) . 

x:dist (x,oll)<3h 

From this it follows that there is a cut-off function a = on f2 having 
the following properties: 

1. < a < 1 on Q. 

2. suppo- CB_U B + . 

3. \u\ = O(ft^) on supp (1 — a). 

4. j n |V<t| 2 < C, uniformly as h — > 0. 

To construct <7, choose \ £ Co°([0>2)), < \ < 1> that equals one on 
[0,3/2] and set u{x) = 1 - \{\x - (±1,0) [/ft) on B± and a = on the 
handle. Properties 1,2 and 4 are easily checked directly, and property 3 
follows from (|39|). 

For brevity, denote, for Q! C fi, 

Qn'(«) = l|V«||i 2(n ,) + a||«||ia(Dnn') 
so that Aq (A) = (u) . Without loss of generality we may assume 

A B (A_) < A B (A + ). 

First, we show 

(40) A B (A min ) > An(A). 

This is easy: Take an optimal pair (u,D) for As(^4 m in), extend u by zero 
outside -£>_ and define a domain D = DL)D' C fi, |-D| = A, by choosing any 



D' C B°i with = A — A m j n . Since u = on D', one gets (40) by using 
(u,D) as a test pair for Aq(A). 

Next, we show a reverse inequality. Using the properties of a and supp V<r C 
supp (1 — a) we obtain, with || • || denoting the L 2 -norm on B±, ||V(ctu)|| 2 = 
\\aVu + (Va)u|| 2 < (||Vu|| + ||Vcr|| max supp(1 _ ff) u) 2 < \\Vu\\ 2 + O(ft^) and 
therefore Qb ± {°u) < Qb ± { u ) + 0(h p ). Now we can use au as test function 
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for the lowest eigenvalue of —A + axDnB± on B±, and this gives the third 
inequality in 

A n (A)=Q Q (u) > J2QB ± (u)>J2QB ± (™)-0(hP) 

± ± 

± Jb ± 

> T\ B± (D ± ) [ n 2 -0(^) 
± Jb ± 

(41) > A B (A_) + (A B (A + )-A B (A_)) / u 2 -0(^). 

Jb+ 

In the last two inequalities we have used property 3. of a, the optimality of 
Ab(A±), and j n u 2 = 1. 

Since we assume Ab(A + ) > A#(A__), this and inequality ([5]) imply 

(42) A B (A min ) > A n (A) > A B (A^) - 0{h?). 

By strict monotonicity of A^ one easily gets from this A_ < ^4 min +o(l) (h — 
0). 

Next, from D C D + UD_llH and \H\ < Ah we have A < A + + A_+4h, so 
A + - A_ > A- 2A_ -Ah> A- 2A min - o(l), and then A min = max(0, A - 
\B1\) < max(0, A — ir) gives 

(43) A + - A„ > mm(A, 2vr - A) - o(l). 

This shows A + ^ A- for h < ho(A, a) and therefore proves part (a) the 
theorem. 

Now we prove part (b). From (^) we have A + —A^ > cq for some constant 
Co > 0, whenever h < ho(A,a), so strict monotonicity and continuity of A# 
imply 

(44) A B (A + ) - A B (A-) >c 

with c > independent of h. Now from ( |40D and (fll|), and using Ab(A-) > 
A_B(^4 m in) (since A_ > A m j n ) and monotonicity, we conclude (Ab(A + ) — 
A B (A^))j B+ u 2 = 0(hP). This and @ give / B+ u 2 = 0(/^).Since, by 

(39), = O(h^), this L 2 bound implies a pointwise bound for u on B + 

by fl47|). Combined with (|39|), applied on the handle, this gives 

(45) sup u(x) 

Finally, we want to deduce from (^) that D c C B- if A > 7r and /i is 
sufficiently small: Since (u, D) is an optimal pair, we have D = {u < t} for 
some t > 0. Equation (45) shows that we are done if we can show that t > c 
for a constant c > independent of h. 

For r € (0, 1) let B_(r) be the closed ball of radius r concentric with B-. 
Applying Lemma ^ to G = 2?_ we see, since ||w||i2^_) > 1 - Oik 13 ) by 
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(ID, that 

(46) inf u > Cr 

B_(r) 

for any r G (0, 1), with c r > only depending on r, A and a, and this 
implies |{u > c r }\ > |5_(r)|. Therefore, we can conclude t > c r as soon as 
|S_(r)| > \Q\ — A. Since < 2tt + 4/i and A > tt, one can find such an r 
ii h < ho, both r and /io only depending on A (and a). This completes the 
proof of the theorem. □ 



4. Free boundary and convex domains 



Proof of Theorem Part (a). First recall, as a consequence of results by 
Brascamp-Lieb [BL] and Caffarelli-Spruck [CS], that the first eigenfunction 
ip on a convex domain possesses only one point where V^> = 0. This point is 
necessarily the point where ip attains its maximum. Now given A, we select 
tn as in Theorem ||[ and we select 5q < tn such that tn + 5o < M where 
M = maxsj rp. With this choice of 5q we use Theorem || to determine a value 
ai for which [Q] tn ~ 5 ° C D C [U] tn+S ° for all a < a\. Then the free boundary 
{u = t} is contained in the closed annulus A = {t^ — 5o < if) < tn + 5$}. We 
have Vip 7^ on A, so C := min^ |V^| is positive. Thus decreasing ct\ to 
a smaller value «o > 0, we can use Proposition [ll| to conclude that for all 
a < otQ we have |Vu| > C/2 on A and hence on the free boundary {u = t}. 
Applying Theorem || we now get the first part of Theorem ^. 

Proof of Theorem || part (b). We only sketch the proof. Fix xq with 
\7tp(xo) 7^ 0. Choose coordinates in which S7ift(xo) = (0,... , 0, a), a > 0, 
and for x' near x' (where x' = [x\, . . . , x n _i)) and t near to = *jj(xo) denote 
the locally unique solution x n of the equation ip(x',x n ) = t by Fo(x', t). For 
a near zero and x near xq one has du a /dx n ^ by Proposition |ll|, so we 
may define .F a similarly for u a instead of ip. 

By a result of Korevaar and Lewis [KL] the level set of ip through xq is 
strictly convex, in the sense that the matrix ( ^..q° . )ij=i,.„ ,n— 1 is positive 
definite at (xQ,io)- Therefore, the result follows if one can show continuity 
of q x in a and (x',t). Now the equation for u gives for a uniformly 
elliptic, quasi-linear equation (writing y = (x',t)) 

n Q^F (v) 

^2 bij(VF a ) a = axG a {yn)Vn - A{a,A)y n 

i,j=i y % ^ 

with bij real analytic and G a = (— oo, t a ], where t a is such that \{u a < 
t a }\ = A. From this it is easy to derive the desired regularity, cf. the proof 
of Lemma 3 in [CGK]. □ 
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5. Numerical results 

In this section we make a few remarks on our method for the numerical 
solution of our eigenvalue problem. 

We use the finite element method for the discretization of our eigenvalue 
problem, with conforming P-l elements. To create the mesh we have uti- 
lized the automatically spatial meshing program encoded by Y. Tsukuda 
(sec [TK]). In order to calculate the approximate first eigenvalue and the 
corresponding eigenfunction, we employ the power method. 

Our method to obtain an optimal configuration is based on an algorithm 
that was introduced in [Pi]. However, we do not insist on D (the sought-for 
optimal configuration) to be a union of elements. This flexibility allows us 
to find a good approximation even without remeshing. 

We now describe the main procedure. The given data are A and a. We 
first take any initial domain Dq satisfying |Dq| = A. Next, if we have ob- 
tained .D n _i (n = 1, 2, 3, • • • ) then we calculate the first eigenvalue A n _i and 
the corresponding eigenfunction it n _i for the finite element approximation 
problem for the operator —A + q%d ii _ 1 . Then we obtain D n from u n _i by 
finding a number to such that < to}| = A and setting 

D n = < t }- 

The number to is determined by a bisection method, i.e. by setting down® = 
0, upo = max(j{4„_i, j = and then iterating Steps 1 and 2 (with L(t) := 
< t}|) 

Step 1: Let interrrij: = (upj + down,j)/2 and calculate L(interrrij). 

Step 2: If L(interrrij) < A, then upj + \ := upj and dowrij+i := interrrij, else 

if L(interrrij) > A, then upj+i := interrrij and downj+i := downj. 

Increase j by one. 

The iteration is stopped when L{interm^) nearly equals A and up^ and 
down/, nearly equal interm^ according to the adopted precision of approxi- 
mation, and then we set to = interrrij. 

Having obtained D n we repeat the procedure above to find u n ,D n+ i etc. 
It is easily seen that A n < A n _i. We iterate until |A n — A n _i| < e, where e 
is given. In the numerical experiments that we have done, we have taken e 
between 10~ 7 and 1CT 10 . 

By the monotonicity of {A n }, the limit limj^oo = Aqo exists. However, 
it is not clear a priori whether Aqo = Afi(a,A) or not. In order to avoid 
the latter case, we have repeated the same procedure with several different 
initial shapes Dq. 

The results of some of the computations that we have done are shown in 
Figures 1-3. They illustrate well Theorems ||, § ||, 0, and |9[ 

6. Some open problems and conjectures 
In this section D = D a ^ will always denote an optimal configuration. 
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Conjecture 1. (Uniqueness and convexity) IfVt is convex then D is unique, 
and D c is convex (at least for a < oq(A) ). 

Concerning the restriction on a compare the remark after Theorem ||. We 
have proved convexity for small a in Theorem ^. 

Problem 1. (Regularity of the free boundary) When is the boundary of an 
optimal configuration smooth everywhere? In general, how can we control 
the size of singular sets of the free boundary ? 

In the convex case we have proved smoothness for small a in Theorem ^. 
A similar method should easily yield smoothness of the free boundary for 
small A and smooth d£l. 

Conjecture 2. In dimension two the free boundary dD is smooth outside 
a finite set. 

We prove some restrictions on the singular set of dD in [CGK]. 

Problem 2. (Topology of D and D c ) If Q is simply connected, is D also 
connected even in the case a > an (A) ( cf. Theorem \^)? If A or \Q\ — A is 
small enough (with a fixed), is D c always connected? 

Compare Proposition |l2| for the case of A close to |f2|. In a dumbbell ipn 
has two maxima. But numerical evidence suggests the following conjecture: 

Conjecture 3. (One component of D c for dumbbell) Let Qh be a dumbbell. 
Then for every a > there is po(a,h) > such that D c consists of one 
component (near one of the maxima of ipn) whenever \Q\ — A < po(a, h). 

Clearly, one would expect po(a, h) — > as a — > 0. 

We now turn to questions of symmetry. A very general problem is the 
following: 

Problem 3. (Symmetry and symmetry breaking) Determine (at least qual- 
itatively) the region in the space of parameters where symmetry breaking 
occurs. 

For annuli the parameters are a, 5 = A/\Vt\ and the ratio r of outer and 
inner radius ('thickness'). For dumbbells the parameters are a, A and the 
thickness of the handle. 

First results on this general problem are given by Theorems || and 0. 
The next three conjectures address other aspects of this problem, i.e. they 
concern other regions in parameter space. They are motivated by numerical 
experiments. 

Conjecture 4. (Symmetry on dumbbells) Let Vt^ be a dumbbell. Then for 
every a > there is p\(a, h) > such that symmetry breaking occurs if and 
only if — A < pi(a, h). 

Conjecture 5. (Symmetry on annuli) For each a, 5 > there is To(a,8) 
such that symmetry breaking occurs for the annulus of thickness r if and 
only if t < To(a,5). 
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Theorem || gives one half of this. The other half means that the optimal 
configuration is rotationally symmetric for 'thick' annuli. Some aspects of 
this conjecture are discussed in [CGK]. 

More generally, it would be interesting to prove symmetry preservation 
in any situation not covered by Theorem Q (i.e. in a non-convex situation). 
In particular, a natural conjecture is: 

Conjecture 6. (Symmetry preservation for small a) For any domain Q 
and any A there is oto(A, f2) such that for a < cto(A, f2) any optimal config- 
uration D has the same symmetries as £1. 

Also, the analysis of the transition between the symmetric and asym- 
metric situations would be interesting, as well as the shape of asymmetric 
solutions for the annulus. 

Problem 4. (Relation between D and the curvature of dVL) Prove that D 
is fat near points where dSl has large positive curvature. 

For example see Figure 1. For a = and A near zero this should be 
not too hard. See [Kl] for the case a = under additional geometric 
assumptions. From this one should obtain the result at least for small a 
and A by perturbation. In [CGK], Thm. 9, we prove in a model case that D 
is thin near a portion of the boundary which has large negative curvature. 

Problem 5. (Limit a — > oo) Consider the restricted minimization problem, 
allowing only such sets D for which D c is a ball. How does this relate to the 
limit a — > oo in our problem? Where does the center of an optimal ball lie? 

This is motivated as follows: Formally, for a = oo the eigenvalue An (a, D) 
equals the first Dirichlet eigenvalue of D c . (The convergence to this value 
as a — > oo is proved in [HH] and [DKM], for example.) Now by the Faber- 
Krahn inequality (see [Ch], for example), the first Dirichlet eigenvalue of a 
domain of prescribed area is minimal if the domain is a ball. So the optimal 
configuration for a large should be close to a ball, at least when A is close 
enough to |0| (so that a ball of volume — A fits into $7). 

Problem 6. (Other Elliptic Operators) Consider the same optimization 
problem for a magnetic Schrodinger operator (iV— o.xdA(x)) 2 with constant 
magnetic field or a uniformly elliptic operator of divergence type — V{(1 + 
oXd(x))V}. 

We have no results for these operators, even if 17 is a ball. 

Appendix: Basic PDE facts 

Here we collect some well-known facts about uniform estimates for solu- 
tions of elliptic equations. We will state these for an equation 

n d 

Pu = 0, P = A + ^2bj(x)- hc(x), x G G, 
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where P has measurable, uniformly bounded coefficients, u € C 1 (G)nC° (G) , 
and G C M n is a bounded open set. In the following estimates, saying that 
the constants depend on P will mean that they depend on sup G (6i, ... ,b n ,c) 
and stay bounded when this quantity stays bounded. 

First, we have the uniform bound (see [GT, Thm. 8.15 and (8.38)]) 

(47) sup |u| < C G ,p{\\u\\ L 2, G ) + sup |«|). 

G dG 
Second, we have Harnack's inequality: If u > on G and G' is a compact 
subset of G then 

(48) sup u < Cq g' p hif u. 

G' G' 

Combining these two we get the following slightly less standard estimate. 
For e > let G t = {x G G : dist (x, dG) > e}. 

Lemma 16. For any e > there is a positive constant cq.p^ such that for 
any u £ C 1 (G) n C°(G) that solves Pu = and satisfies u > one has 

inf-u > cgp M\u\\l*(G) ~ supn). 
G e dG 

Here we set inf u := oo. 



Proof. We have 

u 2 = u 2 + u 2 < \G e \ supn 2 + \G \ G e \ supu 2 

Jg JG e JG\G e G £ G 

< C GtP . e Mu 2 + \G\G e \C GP ( [ u 2 + supu 2 ) 
' ' G t ' JG dG 

where we used Harnack's inequality and the uniform estimate ([47|). If e is 
so small that \G \ G e \ C' G p < 1/2 then we can subtract the last two terms, 
and the claim follows easily. The claim for larger e then follows from the 
fact that infc e , u > inf Ge u if e' > e. □ 
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Figure 1: Ellipse (with a = 1, A = 1, = 6.3). Optimal configuration D 
black. 




(a) a = 1, S = 0.64, r = 3.5 (b) Thinner annulus (r = 1.2) 




(c) Stronger non-linearity (a = 10) (d) Larger area of D (5 = 0.83) 

Figure 2: Symmetry breaking on annuli: The parameters are a (the 'strength' 
of the non-linearity), 5 = A/\£l\ (the relative size of D), and r = r out /r[ n (the 
ratio of outer and inner radius). In each of (b), (c), (d) only one parameter is 
changed, compared with (a). Optimal configuration D in black. 





Figure 3: Symmetry breaking on dumbbells: The parameters are a = 0.1, A = 
1. Symmetry breaking occurs when the width of the handle (2h) is decreased. 
The 'lobes' are unit circles with centers 4 units apart. Optimal configuration D 
in black. 



